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Abstract: Consider a charged Bose gas without self-interactions, confined in a 
three dimensional cubic box of side L > 1 and subjected to a constant magnetic 
field B ^ 0. If the bulk density of particles p and the temperature T are fixed, 
then define the canonical magnetization as the partial derivative with respect to 
B of the reduced free energy. Our main result is that it admits thermodynamic 
limit for all strictly positive p, T and B. It is also proven that the canonical and 
grand canonical magnetizations (the last one at fixed average density) are equal 
up to the surface order corrections. 

1. Introduction 

Much work has been done on the thermodynamic behavior of large systems com- 
posed from independent quantum particles in the presence of external magnetic 
fields. As is well known, the fundamental problem consists in proving the exis- 
tence of the thermodynamic limit for the potentials and the equations of state 
defined at finite volume. In the particular case of the canonical magnetization 
(defined as the partial derivative with respect to the magnetic field of the reduced 
free energy), one has to prove that the derivative (performed at finite volume) 
commutes with the thermodynamic limit of the reduced free energy. Although 
the quantum canonical ensemble is (from the physical point of view) the most 
important one, most of the previous works were carried out either using the 
Maxwell-Boltzmann statistics or in the framework of the quantum grand canon- 
ical ensemble, because in those settings, many physically relevant quantities can 
be expressed employing the integral kernel of the Gibbs semigroup associated to 
the one particle problem. Moreover, one is able to go beyond the bulk terms and 
investigate finite size effects. 

Take for example the grand canonical pressure of a quantum gas in a con- 
stant magnetic field. The rigorous proof of its thermodynamic limit goes back at 
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least to Angelescu and Corciovei [A, A-C]; its surface correction (in the regime 
in which the fugacity is less than one) was obtained by Kunz [K]. As for the 
Maxwell-Boltzmann magnetization, nice results were obtained by Maoris et al 
[M-M-P 1,2]; they wrote down even the corner corrections. Notice that in these 
papers the domain A was allowed to be more general, typically convex with 
piecewise smooth boundary. 

Another result concerning the thermodynamic limit and the surface correc- 
tions for the magnetization and susceptibility of a Fermi gas at zero magnetic 
field was obtained by Angelescu et al [A-B-N 2]. Because this paper motivated 
our work, we are giving some more details about it. 

Firstly, as in our setting, their domain was a rectangular parallelepiped and 
the magnetic field oriented after the third direction. They defined the grand 
canonical magnetization mji(j3,z) (susceptibility xa{P,z)) as the first (second) 
derivative with respect to the magnetic field of the grand canonical pressure at 
B = 0, for all z £ C \ (— oo,— 1]. Their main result can be roughly stated as 
follows: 

i. rriA(/3, z) — 0, Vz £ C \ (— oo, — 1]; 

ii. There exists Xoo(P, z) analytic in C \ (— oo, —1] such that for any compact 
K C C \ (—oo, —1] one has: 

lim sup \xa(J3,z) - Xoo(0,z)\ = 0. 

A->oc zeK 

More than that, they gave even the surface correction for susceptibility and 
proved that this expansion is uniform on compacts. Because the relation between 
the fugacity and the grand canonical average density of Fermi particles can be 
always inverted, they were able to express the grand canonical susceptibility in 
terms of the canonical parameters p and (3. Let us stress that B = and A 
a rectangular parallelepiped were crucial ingredients in [A-B-N 2], the uniform 
convergence on compacts being obtained via a substantial use of the explicit 
formula of the integral kernel of the Gibbs semigroup associated to the Dirichlet 
Laplacian. 

In this paper, we are studying the "true" canonical problem for a Bose gas 
at nonzero magnetic field Bq > (in order to avoid the Bose condensation). 
Using a standard procedure (see [K-U-Z, H]) of deriving the canonical partition 
function from the grand canonical pressure (see (|2.27Jl ). we are able to transform 
the uniform convergence on compacts of the grand canonical magnetization (see 
Lemma^) into a pointwise convergence (/?, p fixed and L — > oo) of the canonical 
magnetization; this result is given in Theorem [2J Moreover, we obtain that the 
canonical magnetization rriL (see (I2.29|H and the grand canonical magnetization 
at fixed average density (see (|2.30jl ) are equal up to the surface order corrections. 

Two natural questions arise: what about Fermi statistics and what about 
higher derivatives with respect to B (the susceptibility for example)? Partial 
answers and a few open problems are outlined at the end of the proofs. 

2. Preliminaries and the Results 

Let A = {x £ K 3 | - -j < Xj < j £ {1, 2, 3}}, L > 1, be a cubic box with its 
side equal to L. Then the "one particle" Hilbert space is TCi.l '■= ^ 2 (A); denote 
with H n ,L the proper subspace of ®" =1 ?ii,i — L 2 (yl ra ) which contains all totally 
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symmetric functions. Denote with Hq.l = C the space with no particles; then the 
Fock space is defined as Tl '■= ©„>o^™-£- O ne can introduce the "number of 
particles" operator Nl as the unique self-adjoint extension of the multiplication 
with n on each TL n .L- 

Assume that the particles (each having an electric charge e) are subjected to 
a constant magnetic field B = Be 3 , which corresponds to a magnetic vector po- 
tential Ba = -j e 3 Ax. If c stands for the speed of light, define ui :— (e/c)B. Then 
the "one particle" Hamiltonian (denoted with H\^{uj)) will be the Friederichs 
extension of the symmetric and positive operator | (— iV-wa) 2 defined on 
Cg°(A). 

Due to the regularity of A, H\^{oj) is essentially self-adjoint on 

D = {/ G C 2 {A) n C^A), f\ dA = 0, Af G L 2 (yl)} . 
The Hamiltonian which describes n particles reads as: 

H ntL {u) = H hL {u)®---®I+--- + I®---®H hL {w). (2.1) 

N ^ ' 

"n" terms 

The second quantized Hamiltonian H l (uj) is defined as the unique self-adjoint 
operator on Tl whose restrictions to H nj L coincide with H n ^{w). 

If T > stands for the temperature and \i e M for the chemical potential, 
then define (3 — k ^ T > and z — cxp (/? /i) (the fugacity), where fee is the 
Boltzmann constant. When working in the canonical ensemble, one considers 
that the bulk density of particles p is constant, therefore the number of particles 
is defined as N{L) := pL 3 . 

As is well known (see [R-S 4]), Hi^l(u) is positive, unbounded and has com- 
pact resolvent; these imply that its spectrum is purely discrete with accumulation 
point at infinity. Moreover, from the min-max principle it follows: 

w£*(H ltL (w)) > wia(H li00 (w)) = |. (2.2) 

It is also known that the semigroup Wl(/3,oj) :— exp {—(3H\^l(u)) is trace class 
and admits an integral kernel G Wi l(x, x'; /3), which is continuous in both its 
"spatial" variables. The diamagnetic inequality at finite volume (see [B-H-L]) 
reads as: 

|G., L (x,x';/3)| < G 0;i (x,x';/3) < ^-1^ exp (2.3) 
If Zi(L 2 (/l)) denotes the Banach space of trace class operators, it follows that: 

\\w L {MU=toW L {M< (2 ^ 3 )3/2 - ( 2 - 4 ) 

Denote with {Ej{w)}j^ the set of the eigenvalues of Hi^l{u)- If /i < 0, the 
grand canonical partition function reads as: 

oo 

S L {0, z,u) = tr^ L exp [-0(H L (u) - fiN L )] = - zexp (-^-(w))]" 1 . 

j=o 

(2.5) 
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The canonical partition function of our system is: 

Z L (/3, p, u) = trn N{L) , L exp (-[3H N{L] ^ L (u))). (2.6) 
The link between them is contained in the following equality: 

oo 

3 L (0,z,u) = ^z n tr Hn>L vq>{-PH n , A ]. (2.7) 
ra=0 

Throughout the entire paper, by log z we shall understand the logarithm function 
restricted to C \ (— oo, 0]. 

Let C be a contour which surrounds the origin, does not intersect the cut 
[1, oo) but contains the spectrum of the trace class operator zWl, where z € 
C \ [exp (/3w/2), oo). Let q(£) = | log(l — £) be an analytic function in the 
interior of C. Define the following bounded operator: 

Q(zW L ) = ±Jd£q(£)(£-zW L )- 1 . (2.8) 

It is easy to see that log(l — zWl) = zWl ■ q(zWL) and using l|2.5(l one 
obtains: 

logSiOS, z, w) = -tr (zW L ■ q(zW L )) . (2.9) 

Employing the above expression, one can easily prove that the grand canonical 
potential (seen as a function of z) is analytic in C\[exp (Poj/2), oo). When \z\ < 1, 
(|2.9I) becomes: 

oo n oo n / p \ 

logS L (/3,*,w) = £— tr(W£) = X)— ( / dxG w , L (x,x;n/3) . (2.10) 

n=l U n=l 71 ^ J A ' 

The grand canonical pressure and density are defined as: 
P L {P,z,w) := ^ log E L (J3,z,w) = -^E^fl-^M) , (2.11) 



and 



61FV 

p L ((3 >Z> cj):=(3z—^(p,z,u>). (2.12) 
Let us remark that pi,(f3,x,Lo) is an increasing function if < x < exp (/3w/2): 

^(/3, x, u) = 1/L 3 tr[(l - xW^)" 2 ^] > 0. (2.13) 
ox 

The proof of the thermodynamic limit for these two quantities goes back at least 
to Angelescu and Corciovei [A, A-C]. Because this result plays an important role 
in our work, we shall reproduce it here. In order to do that, let us define (w > 0): 

i °° 

pm*,») ■■= "Tt^m 5> /2 (^- (fe+1/2) ^) , (2.i4) 
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and 



n n 1 00 

Poo {p, z, uj) := Pz-Q^(fl, *,<*) = P u (2n0 )3/2 E 9i/2 (z e - (fe+1/2)w/3 ) , (2.15) 



(2*/J)' t _ n 
where <?o-(C) are t ne usual Bose functions: 



J> Jo 1 - Ce- 



^(0 = T^TZT / * 1 — ■ ( 2 - 16 ) 



analytic in C \ [l,oo) and if |£| < 1, they are given by the following expansion: 



n=l 

Then the following result is true (see [A, A-C]): 

Theorem 1. Let K C C \ [exp (J3u/2), 00) fee a compact set. Then the grand 
canonical pressure and density admit the thermodynamic limit i.e.: 

lim sap\P L {fi,z,u)-P oo {l3,z,u)\=0, (2.17) 

and 

lim s\xp\p L ([3,z,uj) - p oo (/3,z,uj)\ = 0. (2-18) 
L^oo ^ eA - 

Firstly, because Pt, and are analytic functions, then via the Cauchy inte- 
gral formula it follows that all their complex derivatives admit a limit which is 
uniform on compacts. In particular: 



lim sup 



-5— {P, z, uj) - -^—{p, z, uj) 
oz oz 



0. (2.19) 



It can be seen from l|2.15f) that ]mx x y e p u ,/2 Poo(P, x, uj) — 00, which means 
that the Bose condensation is absent when a nonzero magnetic field is present. 
A very important consequence of the theorem is that the relation between the 
fugacity and density can be inverted for all temperatures and moreover, if < 
Xoo (0, p, uj) < e^l" 1 is the unique real and positive solution of the equation 
poo(P,x,uj) — p and if Xl(/3, p, u>) is the unique real and positive solution which 
solves pl(j3,x,w) = p, then Hhil^ocXl = £oo- 

Let us perform the Legendre transform at finite volume: 

f L (J3, p, uj) := -P L (/3, x L (/J, p, uj), uj) + | log x L {p, p, uj). (2.20) 

A straightforward result is that has the following limit: 

foo{P,p,uj) := -P 00 (P,x 00 (f3,p,uj),uj) + - log Xoo (^p,^). (2.21) 
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Let us denote with ^p{P,^o) the following integral which makes sense in 
the norm topology of B(L 2 (A)) (see [A-B-N 1]): 

drWLiP-r,^) [a- (p-uob)]W l (t,cjq). (2.22) 

Jo 

A particular case of the problem treated in [A-B-N 1] is that wo) is 



even trace class and moreover, for Su sufficiently small one has: 

dW L 



W L (J3, u + 6lu) - W L (fi, uj ) - Suj- 



dui 



(2.23) 



Using for Pl(/3, Z, u>o) the following expression 

1 



3 trpoga-zW^OS.wo))], 



then the estimate (|2.23|) justifies the definition of the grand canonical magneti- 
zation: 



r L ((3,z,Lu ) 



dPr 



(j3,z,u ) 



tr 



(1 - zW L (p,vo)) 



du> 



c duj vr ''~'~ u/ c/3L 3 

(2.24) 

From its definition, one can easily see that Z/, has the same domain of an- 
alyticity in z. Now let us define the natural candidate for its thermodynamic 
limit: 

p BP, 

r^z^o)- w,,). (2.25) 

C OUJ 

Our main technical result is presented in the following lemma: 

Lemma 1. Let K C C\[exp (/3cj/2), oo) be a compact set. Then the grand canon- 
ical magnetization admits the thermodynamic limit i.e.: 



lim sup \r L (P,z,uj) - Too (A z,oj)\ = 0. 



(2.26) 



Let us go back to the canonical ensemble. From l|2.7|l . (|2.11|) and (|2.6|) . one 

can write down an useful representation of the canonical partition function: 



2m 



Ci s 



exp(fP L (/3,£,w) 
I 



N(L) 



(2.27) 



where C\ is a contour which surrounds the origin and avoids the cut. The reduced 
free energy reads as: 



f L (/3,p,uj) := - j^j\ogZ L ((3,p,uj) 
The canonical magnetization is defined as follows: 

m L ((3,p,uj) ■= -^(ftftw). 

C OUJ 



(2.28) 



(2.29) 
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We expect that uil should be close to the following quantity: 

rh L ((3,p,uj) := -^-((3,p,uj) = r L (/3, x L (0, p, w), (J), (2.30) 
c out 

which converges to .Too(/3, Xoo(P, P, oj), u>). 
We are able now to give our main result: 

Theorem 2. Fix < 5 < 1/2. For all strictly positive temperatures, bulk den- 
sities and magnetic fields, the canonical magnetization mj,(/3, p,u) admits the 
thermodynamic limit. Moreover, there exist two positive constants Cg((3, p,lu) 
and Ls(P, p,u>) such that for all L > L$ one has: 

\m L -m L \<C s L~^ 2+s . (2.31) 

Remark. It is clear that mi is a much more convenient quantity. If (at least for 
dilute gases) one would be able to write down an expansion for rfiL which takes 
into account the surface corrections: 

m L — moo + Y m s + °( 1 /^)i (2.32) 

then the estimate l|2.3H would imply that the same expansion is true for mi, 
too. 



3. The Proof of Theorem El 

At this point, we shall consider that Lemma ^ is true and give its proof in the 
next section. In order to simplify the notations, we shall drop the dependence 
on /3, p and lj but we shall reintroduce it when needed. 

The main idea of the proof consists in isolating the principal part of the 
integral from l|2.27|l . Although this procedure is far from being new (in the 
physical literature it is known as the Darwin-Fowler method; see [K-U-Z, H] and 
references therein), we decided to give a rather detailed proof in order to have 
a clearer image of the remainder from Ij2.31|l . Firstly, let us choose the contour 
C\ as follows: 

d :={x L e**, (f>& [-7r,7r]}. (3.1) 
Using 12.20f) . the formula (|2.27() can be rewritten as: 



Z L = exp(-/?/ L L 3 )^ J dcj) exp 



N(L)f3(P L (x L e^)-P L (x L )) 



e -iN(L)4> 



= ex P (~pf L L 3 )N(L)- 1 / 2 A L , (3.2) 
where Al{P, p, oj) is given by: 



Al = V N ( L ) d( p e [N(L)p/p(KP L (x L e'*)-P L (x L ))] e iN(L)[0QP L (x L e'<l>)/p-<t>]_ 

(3.3) 

Denote with p o(</ ) ) the function given by 5RP 00 (x 00 e^), and with pl(4>) the 
function given by 3?Pf,(a;i,e*^). 
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,2 ~ 

We shall prove that /j? (0) < and that the following limit is true: 



lim Al ■= Aoo = . /-- 



L — >oc 



27T/3 



(0) 



(3.4) 



Firstly, let us remark that Pl(£) = Pl(Oj where the over line means complex 

0u>/2 _ 

conjugation. Let < e < - 2 x °° be a fixed positive number, and let C e be a 

circle centered in Xoo with radius e. If £ belongs to the interior of C e then the 
Cauchy integral formula gives (L < oo): 



(3.5) 



Take < 5 < 1/2 and define <f> L := A^L)" 1 / 2 ^/ 3 . Then for L big enough and 
\<t>\ < <Alj one nas that xlz 1 ^ belongs to the interior of C e /2 and l|3.5|l implies: 



pl{4>) 



Atti 



d(P L (0 



(3.6) 



Theorem ^implies now that: 



Hm^(0) = ^(0) and V|0| < 0; 



< const. 



(3.7) 



Denote with Sl(4>) the function SP^^Le 1 ^). In a similar way, one can obtain 
that for all [(f>[ < </>l: 



"(0 



< const. 



(3.8) 



Finally, let us prove l|3.4|) . This will be done in two steps: 



1. The first one consists in showing that the contribution to Al coming from 
the region \<fi\ > <t>L is exponentially small in L; we shall use that pl{4>) is an 
even function, is decreasing on the interval [0, tt] and has a non-degenerate 
maximum in 0; 

2. The second one consists in a more careful study of the case in which \(f>\ < 4>l\ 
the most important thing here is showing that the oscillations of the imaginary 
part are small. 



Let us prove 1. Firstly, because Pl(4>) = Pl{- 
at <fi £ [(f> l, 7r]. From (|2.11|) one has: 



it is sufficient to look only 



Pl{4>) 



1 



2f3L 3 



^log (l-a; L cos0e f)Ej ) 2 + x\ sin 2 



-2{iE 3 



(3.9) 



and: 



dp L , ,s xl sin 



f3L 



'• ^ |1 -x L e x <l> e-^| 2 



(3.10) 
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This means that pL is decreasing on [0, 7r]. Its second derivative in reads as: 
d2 P L (Q) - XL V e ~ PEi mi) 

For L big enough, one has that xl > £oc/2. The above equality implies: 

f^ (0 ) < -||>W L . ,3,2, 
To obtain an uniform estimate in L, we employ (see [K]): 

J-trW L = 1 3/2 ■ Hfa/o\ 0- + °( 1 / L ))- ( 3 - 13 ) 

L 3 (2tt/3) 3/2 smh(w/3/2) 

Then IpTTjl . (|3~T!i|) and l|3~T3)l imply that for L sufficiently large: 

-rpr (0) < — J3" and — 77^(°) < °- ( 3 - 14 ) 

d<j) 2 4/3 (2tt/3) 3/ smh(w/3/2) d</> 2 

From (|3.14() and l|3.7|) it follows that for all <f> E [4>L,n] one has: 

p L (0) - P L (x L ) < PL ) ~ PL (0) < \ ^ (0)01, (3.15) 
or in other words, 



exp( N ^ \p L (<l))-P L (x L )] \ <exp r 



(3.16) 



therefore the contribution to Al coming from \4>\ > is exponentially small. 

Now let us study the region \ tp\ < c/>l- We are mainly interested in the behavior 
of the imaginary part of the pressure (see l |3.8| l). Similarly as in i|3.9|) . one has: 

1 » xl sin <j> e~@ Ej 



SL ^ ] = ~JU^ arg (1 " x ^ e ~ PBi ) = JljY. arctan i 



L 3A^ & ^ > (3L 3 ^ l-x L cos</>e-' 3 ^' 

3 3 

(3.17) 

By direct computation, one obtains: 



Together with (|3.8(l . one has: 

/3 -«L(0)-^ = O(^i). (3.19) 



We are interested now in the following integral: 



V N ( L ) f' PL d e {N(L)f3/p(3tP L (x L e'*)-PUxL))] e iN(L)lf3^P L (x L e*<t>)/p-<l,}^ 
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Changing the variable in t = y/ N(L)<fi one obtains: 



1 r N ( L ) s/3 



2tt 



dt e 2 "' 



(0) t 2 e G[AT(L)- 1 / 2 + 5 ] 



JV(L)V3 



(3.21) 



Using (|3.14|) , (|3.7|) and the Lebesgue dominated convergence theorem, it is easy 
to see that the above integral converges to (|3.4() . 

From (|2.28|l and (j3.2(l . the reduced free energy reads as: 



h = f L + 



lQgjVQC) _ log^L 

2/3L 3 /?L 3 



(3.22) 



From (|^2"^|) and l(2~3U|> it follows: 



e 1 &4 L 



(3cL 3 A L duj ' 



(3.23) 



m L - m L 
_ N(L) 1 / 2 
~ 2tt A l ., , 

where: 



e d 
c dio 



(3.24) 



(P L (x L e 1 *) - P L (x L )) 



--^ (PUxLe^) - P L (x L )) 



r L (x L e*) - r L (x L ) - - 

c 



dz > duj 



(3.25) 



dP L . . dx L 



From $TT% it follows that - £ ^ (/3, 2, cj) = f3z^((3, z, oj). This implies 



C OLU OZ 



dp L , , 



Lemma ^ and Theorem ^ imply (via the Cauchy integral formula) that the 
above quantity remains bounded when L goes to infinity. Performing a similar 
analysis of the integral from H3.24J) as that one made for Al, it follows that 
the contribution coming from the region \<f>\ > 4>l is exponentially small in L. 
When we are analyzing the region \<j>\ < </>£, the factor N(L) 1 / 2 disappears when 
we are changing the variable, therefore wil — rhh will behave as the quantity 
from (|3.25|l . and this one behaves like 4>l ~ L~ 3 / 2+s . The proof of (|2.31|l is now 
completed. 
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4. The Proof of Lemma ^ 

The idea of the proof is borrowed from [A]. Namely, the uniform convergence 
on the compacts which belong to T> := C\ [e' 3 " / 2 , oo) will be obtained applying 
the Vitali Theorem to the sequence of analytical functions Zz,(/3, z, ujq) when L 
is going to infinity. Therefore, we have to make the following three steps: 

I. Prove the uniform boundedness in L of the functions /~z,(/3, z, u)q) on any 
compact K C T>; 

II. Identify the limit (-Too(/3, z, u>q) in our case) and prove that it has the same 
domain of analyticity T>; 

III. Prove the existence of a set T>q C T> having at least one point of accu- 
mulation, such that .Tl(/3, z, u>q) has a pointwise convergence to /^(/J, z, ujq) on 
■D . 



4-.1. The proof of I. Here is the "nontrivial part" of our paper. In order to see 
where the main difficulty is, let us take a look at l|2.24[l . One could try a bound 
on Jl(/3, z, uj) using the following inequality (see also l|2.22(l ): 



sup 



tr 



(1 - zW L ) 



< 



dWr 



<C[fi,K,u) 



sup || 
dW L 



(i - zW L y 



B(L 3 ) 



du> 



Unfortunately, the linear growth of the vector potential a leads to a bad trace 
norm estimate for ^P^- , which behaves like L A and not like L 3 as needed. In the 
case of the pressure [A-C], this kind of problem does not appear (see (|2.9(l . ()2.8|l 
and l|2.4|l ). because for the Gibbs semigroup the trace and the trace norm are 
equal and grow like L 3 . 

Because the proof is rather lengthy, we shall outline in a few words our strat- 
egy. Firstly, having 1)2. 8JI in mind, let us define (wq,/3,t > 0): 



(4.1) 



where ^ belongs to the contour C. Then the pressure will admit the following 
representation (see i[2.9|l ): 



P L (f3,z,u ) = - 



1 

2m 



dtq(t)jptr\g Uo (t,z;p,p)], (4.2) 



where the norm convergent integral from (|2.8|) commutes with the trace from 
(|2.9(l . Let us argue that we can choose the same contour of integration if u> varies 
in a small interval f2 := [u>q, lo\\ 1 < uj$ < uj\. In other words, we search for a C 
included in C \ [1, oo) such that for all L > 1, z € K, u> € Q and £ G C one has 
the following bound: 



m-zWLi^u)]-^ <M<oo. 
Assume that every z G K verifies the following condition: 



(4.3) 



dist{z, [ e 2 ,oo)}>5>0, \z\ < d < oo. 



(4.4) 
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We know that the spectrum of Wl(/3,(j) is included in [0, e"^ 1 ] for all L > 1 
and oj S ft. We claim that C can be chosen as the union C\ U C2, where Ci is 
given by (rj > 0): 

{(l + t,±i7)| - »7 < * < 2d} U {(1 - 77, *) 1 -77 < t < r)}, 

and Ci is chosen such that if £ G Ci, then |£| > 2d + 1. It is not difficult to prove 
that by choosing 77 sufficiently small, then: 

sup sup sup |(£ — z r) _1 | < M < 00, (4.5) 

s 0<r<e~ 2 

and via the Spectral Theorem, (|4.3[l takes place. 

If we manage to prove the existence of a numerical constant c(/3, if, Wo) such 
that for all 6u> S (0, wi — Wo) and 2 € K to have: 

supsup 7 J-|Pl(/3,z,w + 5w) - P L (J3, z, u ) | < c(/3,K,w ), (4.6) 

then the magnetization will be bounded by the same constant. This estimate is 
straightforward if a stronger one takes place: 

supsupsup — ^— \tr[g Uo+gu (£,z;0,/3)] - tr[ ffwo (£,z;/3,/3)]| < C(p,K,u) Q ). (4.7) 



Our main task will consist in constructing a trace class operator -A^o-H^CC) z; (3) 
having the following two properties: tiAu +Su(£, z ] 0) = tr<j Wo (£, z; /3, /?) (i.e. its 
trace is not depending on Sco) and moreover, 

supsupsup— ^\\g uo +6w(€,z;/3, (3) - Aw +8u{€, z; 0)\\ii < O'((3,K,u} ), (4.8) 



z £ 5u 



which would clearly end the problem. 

We will see that (u = ujq + 8lu) A u (£, z;0) can be chosen as a product 
£fw(£ 5 z ', (3, (3/2)S UJ {f3/2) where the first term is bounded, the second one is trace 
class and moreover: 

sup sup \\g u (£, z- (3, /3/2) - , z; (3, /3/2)|| B(L 2 ) < Ci(fi, K, cJ Q )Scj, (4.9) 
and 

\\W L (fi/2,u) - S u (0/2)\\ Xl < C 2 (/3,^ )i 3 fe. (4.10) 
In particular, l|4.9|l and (|4.10|> imply: 

\\S u (fi/2)\\ Xl < C(I3,lu )L 3 and sup sup \\g u (C, z;/3,0/2)\\ fl(LS) < C 3 (p,K,wo). 

(4.11) 

Employing these estimates together with the following identity (see (14. If) ): 

g u (t, z- (3, (3) = (£, z; /3, /3/2)Wj (/3/2, w), 

the proof of l|4.8|l follows easily. 

The rest of this subsection is dedicated to the rigorous proofs of these esti- 
mates and will be structured in a sequence of technical propositions. We start 
with a well known result, given without any other comments: 
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Proposition 1. The Dirichlet Laplacian defined in A admits a trace class semi- 
group which has an integral kernel given by the following formula: 

3 

G ,l(x,x';/3) = l[g 0iL (x j ,x' j ;(3), (4.12) 

where the "one dimensional" kernels read as: 
9o,l{x,x';0) = 

= (2^7* £{ 

:= go, oo (x,x';0) + to, L (x,x';0). (4.13) 



(x-x' + 2mLy 



2/3 



cxp 



(x + x' - 2mL - Ly 



Using the previous proposition, one can write: 

G ,l (x, x'; 0) = G ,oo (x, x'; 0) + Z ,l (x, x'; 0), 

where: 



G ,oo(x,x';/3) 



(2tt/3) 3 / 2 



cxp 



,'|2 



2/3 



(4.14) 



(4.15) 



The purpose of the next proposition is to give a few properties of smoothness 
and localization of the reminder Z 0i l(x, x'; 0): 

Proposition 2. For all > and L > 1, there exist two positive numerical 
constants C\ and c 2 such that: 



dZ, 



0,L 



(x,x';/3) < ci i^y^ G ,oo(x,x';c2/3); 



n. max 



0,L 



80 



d 2 Z ,_ 



dxjdx' k 



(x,x';/3), 

hi. \Z 0iL \(x, x'; /?) < ci (1 + 0) G ,oo(x, x'; c 2 /3). 
Proof. For x,x' e (-L/2.L/2) define: 

Cl(x '"' ;/?) = 7^ S exp ( 



(x,x';/3)| < Ci^pG ,oo(x,x';c 2 /3); 



ez\{o} 



(x - x' + 2mL) s 
2£ 



(4.16) 



C 2 (x,x';/3) = - 7 ^^exp(- 



[x-x'- (2m+l)L] 2 
2^ 



It is clear that if one obtains uniform estimates in L and for these two 
quantities, the same would be true for Zo,io too. A very useful estimate is the 
following: 

Vt > 0, te _ * = 2(i/2)e-*/ 2 e- t / 2 < 2 e - t / 2 . (4.17) 



14 



H. D. Cornean 



Then we have the inequalities: 



dx 

d 2 g 

dxdx' 



, . . const. v 

(x, x;P)< ex P 

P m=£0 



(x - x' + 2mLf 



(x, x'; 13) < ^2 exp f - 



^3/2 



m#0 



4/3 

(x - x' + 2mL) 2 



(4.18) 



4/3 



and similarly for £ 2 - At this point we have to control the summation over m. Let 
us prove the following inequality: 



^exp 



{x-x' + 2mL) 2 \ ( {x-x'f 

1 < ci(l + /3) cxp - 



4/3 



c 2 /3 



(4.19) 



Because \x — x'\ < L one has (\m\, L > 1): 

(x-x' + 2mL) 2 = {x-x') 2 +AmL{x-x')+Am 2 L 2 > {x-x'f : +4(|m|-l). (4.20) 
Therefore: 



ex P 1 - 



(x - x' + 2mLf 



4/3 



< 



< exp 



2 exp 



(x-x'f 
( (x - x') 2 



1+ ^cxp(-im) 

m>l 



! 2/3 



sinh 



(A) 



2/3 exp (--^ 



V 2/3 



< 



< 2(1 + /3)exp 



(x-x') 2 
4/3 



In order to control £2, one has to study the following quantity: 
^ / [x + x'-(2m+l)L] 2s 

A ■= E «p ( ^ 

m£Z 



(4.21) 



(4.22) 



Denote with £ = x + L/2 and £' = x' + L/2: then < £, £' < L, x - x' = i - C 
and (f + C) 2 > K - O 2 - It follows: 

[x + x - (2m + 1)L] 2 = [£ + £'- 2(m + 1)L] 2 (4.23) 
= (£ + C') 2 - 4(m + \)L(i + £') + 4(m + 1) 2 L 2 . 

If m < —1, then: 

[x + x' - (2m + 1)L] 2 > (x - x') 2 + 4(|m| - 1). 

If m = 0, then: 

[x + x'~ (2m + 1)L] 2 = [(L/2 - x) + (L/2 - x')] 2 > (x - x') 2 . 



On the Magnetization of a Charged Bose Gas in the Canonical Ensemble 



15 



If m > 1, then: 

[x + x' - (2m + 1)L] 2 > (£ + C) 2 + 4:L 2 (m + l)(m - 1) > (x - x'f + 4(m - 1), 

and we can repeat the summation procedure used in 14.19[) . Putting all these 
things together, the proof is completed. 

The next proposition is a variant of the perturbation theory for self-adjoint 
Gibbs semigroups (see [H-P]). Instead of starting with a perturbation of its 
generator, we start with an approximation of the semigroup. Although simple, 
this proposition contains the main technical core of our paper. 

Proposition 3. LetJi := ~L 2 {A) and let H be a self-adjoint and positive operator 
having the domain D. Fix (3q > 0. Assume that there exists an application < 
/3 < /3o — S{P) E B{TL) with the following properties: 

A. sup 0</3 <^ \\S{P)\\ < ci < oo; 

B. It is strongly differentiable, RanS'(/3) C D and s — lim^o S([3) = 1; 

C. There exists a normly continuous application < (3 < (3q — > R(f3) G B(TL) 
such that \ \R{P)\\ < c 2 //3 a where < a < 1 and: 

— f + HS((3)f = R(f3)f. (4.24) 

Then the following two statements are true: 

i. The sequence of bounded operators (n > [1/p]): 

p/3-l/n 

T n {fi) := / dr exp [-(/3 - t)H]R(t) 

Jl/n 

converges in norm; let T((3) be its limit; 

ii. The following equality takes place in B{TL): 

e X p(-[3H) = S(0)-T([3). (4.25) 

Proof, i. The norm convergence is assured by the integrability condition imposed 
on the norm of R(f3). Moreover, when (3 is near zero: 

sup||T„(/3)|| ^(a)/? 1 "", (4.26) 

n 

therefore the same thing is true for T(j3). 

ii. Let < (5\ < (i < /?n. If n > [l//3i] and 4> &T~t, define the vector: 

MP) := exp {-(3H)cj> - S{p)4> + T n [fl)<f>. (4.27) 

From (|4.26(l and condition A it follows: 

limV>„(/3) := ip(/3) = exp {-PH)<f> - S(/3)<j) + T(0)(p and sup \\ip n (/3)\\ < const. 

(4.28) 

Define f n {(3) = 1 1 tAt-i 1 1 2 and f((3) = ||^(/3)|| 2 - From the strong convergence to 
one of S(f3) when (3 goes to zero and from the norm convergence to zero of T(/3), 
it follows that lim^o /(/?) = 0. If we manage to prove that /(/?) is decreasing, 
then it would be identically zero and this would end the proof. 
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Notice that T n {j3) is normly differentiable and: 

-^0 + HT n (f3)<f> = exp ( -~H)R((3 - l/n)0. (4.29) 
The positivity of H implies: 

|| < 2 R[<^»(/3),exp UlifW/J - l/n)<f> - R{[3)<t>)]. (4.30) 



one has: 



Using /„G8) = /nC9i) + I^dr^, 
< /„(/?) < /„(/3i) + const ^dr ||exp R(r- l/n)4>- R(r)(t>\ |. (4.31) 

Because the above integrand is bounded on the domain of integration and has 
pointwise convergence to zero, the dominated convergence theorem implies that 
the whole integral converges to zero. Taking the limit, it follows that f(f3) is 
decreasing and therefore is identically zero. 



Define the "magnetic phase" : 

^(x, x') = x • a(x') = ie 3 • (x A x). (4.32) 

We shall use this quantity as a local gauge transformation (see [C-N] for another 
use of this idea); namely, it will alter the magnetic vector potential by making 
it depend on x — x', only: a(x) — a(x — x') = V x <y9(x, x'). To see how this 
transformation acts (at a formal level) on Hi l(lu) : let us notice the following 
equation: 

[-iVx - wa(x)] e ^( x - x ') = e «"¥>(x,x') [_,v x - ua(x - x')] . (4.33) 



Proposition 4. The bounded operator S(f3) € £>(L 2 (j1)) given by the integral 
kernel e luJl P( x < x )Go,l(x, x'; (3), verifies the hypotheses of Proposition^ having 
the following properties: 

i. s - lim/j^o S(j3) = 1; 

ii. The application (0,oo) 3 (3 i — ► S(f3) is strongly differentiable and uni- 
formly bounded; 

iii. RanS(/3) £ Dom(H 1<L (cj)) and §§(/?)/ + H 1<L (u)S(P)f = R{fi)f, where 
R(f3) has an integral kernel i?(x, x' '; (3) given by: 

e"^ [w 2 a 2 (x - x')G ,l(x, x'; 0) + 2«^a(x - x') • V x Z , L (x, x'; 0)}. 

Proof, i. Throughout the whole section, we shall use a few well known results, 
which are given without proof. Let us begin with an useful boundedness criterion 
for integral operators (see [S]): 



On the Magnetization of a Charged Bose Gas in the Canonical Ensemble 



17 



Let A be an integral operator, given by a continuous integral kernel A(x, x') € 
C(A x A). If the next inequality holds: 



sup / dx' |A(x,x')|, sup / dx |A(x,x')| 



< C < oo, 



(4.34) 



then the operator norm of A in B(L 2 (A)) is bounded by C. 

If c, p and r < P are three strictly positive numbers, then we have the 
following three identities: 

3 



dy exp 



|x-y| 2 |y-xf 



c(/3 - r) 



CT 



p 



exp 



,'|2 



c/3 



dy 



■ exp 



|x-y| 

c(3 



= 1, 



dy |G ,oo(x,y;/3)| ; 



(27T/3)" 



A very useful inequality will be the next one (t > 0, n > 1): 



T exp 



/9 



< const(n)/? 2 exp 



2/? 



(4.35) 
(4.36) 

(4.37) 
(4.38) 



Let us get back to the proof of i. Firstly, because \e luMp ( x ' x )Go,l(x, x'; f3)\ < 
Go,oo(x, x'; /3), it follows that S((3) obeys the condition (|4.34() with C < 1, which 
means that is uniformly bounded in [3 > 0. 

Let us show now that the operator given by the integral kernel (e l " v ( x ' x ) — 
1)Go,l(x, x'; pi) converges in norm to zero. Let us remark first that \<p(x, x')| < 
L |x — x'| and moreover: 

|(e^ x ' x ') -1)| <w|^(x,x')| <wL|x-x'|. (4.39) 

Then (using with n = 1): 

|( e -^ x ')_ 1 ) Go L ( XjX ' ;/3 )| < const L/3 1 / 2 G ,oo(x,x'; 2/3), (4.40) 

therefore its operator norm behaves in zero like p 1 ^ 2 at least. The proof of i is 
now straightforward. 

ii. We will prove that the application is in fact normly differcntiable. For 
p > and SP sufficiently small, one has: 

S(P + SP)f-S(p)f = 6P f dx' e ^(- x ')^i(.,x';/3)/(x') + 



(SP) 



2 ' a 2 G„, 



dx' e ^(-' x ')- 



-(.,x';P)f(x% (4.41) 



dp 2 

where p is situated between P and P + Sp. It is not difficult now to see that the 
"operator derivative" is an integral operator whose kernel is the derivative with 
respect to P of the initial one. 



18 



H. D. Cornean 



iii. Denote with Dq the common domain of essentially self-adjointness for 
H hL {u)),u>0: 

D = {^e C 2 (A) n C 1 ^)! i>\ aA = 0, e L 2 (yl)}. (4.42) 

The action of H\^(w) on a function from Dq is as follows: 

[fli iL («)^](x) = -(^)(x) + 2zwa(x) • (W)(x) + w 2 a 2 (x)V(x). (4.43) 

Now take / € Cg°(.A). After integration by parts, using (|4.33(1 and the fact that 
Go,l(x, x'; (3) solves the heat equation in the interior of A, one obtains: 

{H liL (u)il>, S(P)f) = [ dx! dx Wif^V^^ ■ 



• [-Ac + 2zwa(x - x') • V x + w 2 a 2 (x - x')]G ,l(x, x'; (3) 
= -(^,S'(P)f) + (^,R((3)f). (4.44) 

The result follows easily after a density argument and with the remark that: 

a(x-x')- V x G ,oo(x,x';/3) = 0. (4.45) 

Finally let us notice that the norm of R{(3) is independent of L and is inte- 
grable in zero. To do that, one has to employ the estimates from Proposition 
(|4.38j) and the criterion from Ij4.34jl . 

In order to perform a similar perturbative treatment of the semigroup near a 
nonzero magnetic field, we need the estimate given by the next proposition: 

Proposition 5. Let to. be a unit vector in R 3 . Then there exist three positive 
numerical constants s, C4 and C5 such that for all lj £ fl, x, x' £ A and (3 > 0, 
one has the following uniform estimate in L: 



(I + 3Y 

|n ■ (-tV - a;a(x))G M , £ (x, x'; 0)\<c A y ^ exp 



Proof. Proposition |21 allows us to write down the following integral equation 



(4.46) 



W l (J3,lj) = S(/3) - drW L (f3~ t,lu) R(t). (4.47) 
Jo 

Because S(f3) and Wl are self-adjoint, one can rewrite (14.47ft as: 

W l (P,oj) = S(/3) - dr R*(t)W l (P-t,lu). (4.48) 



In terms of integral kernels, (|4.48|l reads as: 

rP 



G w . L (x,x';/3) = S(x,x';/3)- f dr f dy i?*(x, y, T)G u>L (y, x', ; /3-r), (4.49) 

JO J A 

where the equality is between continuous functions in C(A x A) and the integral 
in t has to be understood as "J"' 3 e " in the limit e \ 0. 
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Because the kernel of R* reads as i£*(x, y; r) = i?(y, x; t), by direct compu- 
tation one can obtain the estimate (see Proposition^: 



fl + t) s ' 

|n- (-zV-o;a(x))i?*(x,x';r)| < c' 4 - ^— exp 



(4.50) 



where one has to apply 14.3311 then use the estimates from Proposition^] finally, 
introducing H4.50jl . (|2.3jl . (|4.35(l and (|4.38jl in Ij4.49jl and because the singularity 
in t is integrable, the result for G Wi l is straightforward. 

Take uj = uj + 8lo G fl. The analogous of Proposition 0] at nonzero magnetic 
field is: 

Proposition 6. The bounded operator denoted with S u ((3) and given by the ker- 
nel e ll5 ^( x < x )G W0! l(x, x'; (3) has the following properties: 

i. (0,oo) 3 (3 i — > S u ((3) is strongly differ 'entiable and s — lim^o S u (l3) = 1; 

ii. RanSU/3) G Dom(H 1<L (u)) and $jS u (0)f + H ltL (u)S u ((3)f = RUP)f, 
where R u ([3) is given by: 

iC(x,x';/3) = e <^)^') [(fe) 2 a 2 (x - x')G Wo , L (x, x'; 0) + 
+ 2(fc)a(x - x') • (zV x + ^ a(x))G Wo , L (x, x'; /?)] . (4.51) 

Proof, i. Rewriting the integral kernel of Sa,(f3) as: 

^(x,x';/3) = G w , L (x,x';/?) + ^M^') _ i] G w , L (x, x'; /3), (4.52) 

and using the diamagnetic inequality, one can reproduce the argument from 
(14 .40(1 in order to prove that the second term converges in norm to zero. Clearly, 
the first one converges strongly to one. 

If {ipj} and {Ej} denote the sets of eigenvectors and eigenvalues of Hi,l(w), 
then: 

G^ L (x,x';/3) =^e-^Vi(x)W(x / ) ) (4-53) 

3 

where the series is absolutely and uniformly convergent on A x A. This can be 
seen from the fact that the semigroup is trace class and that the eigenfunctions 
belong to Dq and admit the estimate: 

iV'il(x) < const(L) (Ej + 1), (4.54) 

obtained from the fact that the resolvent [#1,1,(0;) + is bounded between 
L 2 (A) and L°°(yl). It follows that uniformly in A: 

\G„A-,-;P + 8P)-G u , L (-r;P)--^(-,-;P)\ < const(i) (5/3) 2 , (4.55) 

which is sufficient for the strong differentiability (see also l |4.41|l '). 

ii. One has to make the same steps as in the proof of the third point of 
Proposition^] As for the norm of R(/3), let us see that is independent of L and 
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is integrable in zero: from Ij4.51|l . (|4.46[1 . (|4.38(l and (|2.3(1 . one can obtain an 
estimate on the kernel of i? w (0) of the following form: 

x'; 0)\ < cqSu)(1 + /3) s G ,oo(x, x'; c w 0), (4.56) 

which implies that its i?(L 2 (/l)) norm is bounded by a constant multiplied with 
Suj (see fOU). 

We shall give now without proof a result which gives sufficient conditions for 
an operator defined in B(L 2 (A)) to be trace class: 

Proposition 7. Let {T n } a sequence of trace class operators, converging to T 
in B(L 2 (A)). If sup n HT^I^ < c < oo, then Teli and ||T||zi < c. 

Remark. Assume that an operator T is defined by a i?(L 2 (/l))-norm Riemann 
integral on the interval [a, b], with a continuous trace class integrand S(t). If 
/ dt | l'S'(t) 1 1^ < c < oo, then T is trace class and 



\T\\ Xl <c, tiT= [ 

J a 



b 

dttv S(t). 



Denote with R(uj,0) the bounded operator given by the kernel: 

R(x, x'; 0) = 2 a(x - x') • (-»V - w a(x))G^ L (x, x'; 0). (4.57) 
Among other things, the next proposition proves l|4.10[l : 

Proposition 8. Take > and uj = ujq + Suj 6 Q . 

i. The operator S u! {0) is trace class and moreover, there exists a positive 
numerical constant c such that: 

\\W l (0,u)-S u (P)\\bwa))<c6w and (4.58) 

\\W L (P,u) - S u (P)\\ Xl <cSloL 3 . (4.59) 

ii. For all x, x' € A and uniformly in L, one has: 

G u>i (x, x'; 0) = e 4,5 ^ x < x '>G^ L (x, x'; 0) + (4.60) 
+ Suj f dr [ dy e ^( x ^G W0 , L (x,y;/3-r) e ^^ x '' J R(y,x';r)+O((M 2 )- 



JO J A 

Proof, i. We know that as bounded operators: 



W L (0,u J ) = S iU {0)- dTW L (0-T,cj)R w ( T ). (4.61) 
Jo 

We have already seen that the B(L 2 (A)) norm of Ru{t) is bounded by a 
constant multiplied with Suj (see H4.561 ). Its Hilbert-Schmidt norm is bounded 
by (see (|Q7jl : 

H^WIk < const Suj (1 ^ L 3 / 2 . (4.62) 
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For the semigroup, we know that ||Wi||s(L2(vi)) < 1 and from (|2.3(l and 14.37(1 
it follows that: 

\\W l (J3-t,w)\ |x 2 < const _ 1 t)3/4 (4.63) 

With the help of the well known inequality ||j4 B\\% 1 < ||A||x 2 ||-B||x 2 , it follows 
that in both situations (B(L 2 ) and 2i(L 2 )) the singularities in t are integrable 
(see the previous remark), and the desired bounds follow easily. 

ii. The formula (|4.60() is obtained from l|4.61|l by isolating the term which 
contains Slu. 



The next proposition imposes sufficient conditions on a trace class integral 
operator such that its trace to be equal to the integral of the kernel's diagonal 
(see [R-S 1]): 

Proposition 9. Let T E Ti(L 2 (A)), given by the integral kernel T(x, x') <E 
(7(2 X I). Then: 

trT = / dxT(x,x). (4.64) 

J A 

The rest of this subsection is dedicated to the proof of l|4.9|) . Fix 0, r > 0. Let 
< loq < W i and let [2 — [u)q, u>\\. Clearly, g u (£, z; (3, r) is trace class and admits 
a continuous integral kernel given by the following series, which is absolutely 
and uniformly convergent on A x A (see also (I4.53|H : 

T w (x ) x / )=^K-^e X p(-/3E i )]-^e X p(-r^)^'W^(x') ! (4.65) 

3 

Notice that in order to simplify the notations, we did not specify the dependence 
on £, z, P and t. 

Let us start with the equation satisfied by T^(x, x'): 

Proposition 10. As continuous functions: 

^(x,x')-z/ dyG L ,, L (x,y;/3)T w (y,x') = zG„, L (x,x'; r). (4.66) 

J A 

Proof. The above equality is nothing but the rewriting in terms of integral kernels 
of an identity between bounded operators: 

[{-zW L (p,w)]g u (S,z;/3,T) = zW L (T,u). (4.67) 



For further purposes, we shall prove that for all L > 1, |T w (x, x')| ~ e a ' x x > 
for some positive a. We need first a few definitions: let p(x) := (1 + x 2 ) 1 / 2 and 
a > 0. It is known that the partial derivatives up to the second order of p are 
bounded by a numerical constant and moreover: 

Vx € R 3 , e ±"p(*) e T«|x| < const ( a ). (4. 68 ) 

Fix x G 7l. Denote with A(a) the multiplication operator with e~ ap ( '~ x °). Then 
A(a) and A _1 (a) = A(— a) are bounded operators and invariate Dq (see I4.42t ). 
An useful result is contained in the next proposition: 
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Proposition 11. i. The operator A(a)Wi(r, oj)A(— a) belongs to B(L 2 ), and 
has a norm which is uniformly bounded in L, Xq € A, u> € £2 and < r < (3; 

ii. The operator A{q)Wl{!3, uj)A(— a) belongs to i?(L 2 ,L°°), having a norm 
which is uniformly bounded in L, lo and xo. 

Proof, i. The inequality (|4.68l) allows us to replace A with the multiplication 
operator given by e Q ' x_x °'. Let ip € L 2 (/l) and define: 

<f>:=A(a)W L (T,cu)A(-a)= [ dy e - Q '(- x °>G L ,„(-, y; T)e a "^'H{y). 

J A 

(4.69) 

Let us remark an elementary estimate, which is true for all < r < (3: 

< const(a,/3). (4.70) 
Applying (|4.t)8|) . the triangle inequality, (|2.3|) and 14.70f) one obtains: 
|0|(x) <const(a) / dx' e Q l x - x '' \G u<L (x, x'; r)| |^|(x') < 

J A 

<const(a,/3) / dx! G ,oo(x, x'; 2r)|V|(x'). (4.71) 

J A 

The result follows from jOS} and l(Q3|) . 

ii. We apply the Schwartz inequality in l|4.71|l with r = f3 and then use (|4.37|) . 



e «|x-x | exp 



Proposition 12. Under the same conditions as above, there exists a sufficiently 
small < a < 1 such that the following inequalities are true in B(L 2 (A)), 
uniformly in L, to and x : 

i. \\Wl(P,u) - A(a)W L (l3,uj)A(-a)\\ < a const 

ii. Uniformly in £ € C and z G K one has: 

\\A(aM ~ zWzOMr^-tOlkL'M) < const(/3). 

Proof, i. Let S(j3) = A(a)W L (j3,u)A(-a). We will see that S(j3) obeys the 
conditions of Proposition [3] From Proposition ^] follows condition A. Then 
S([3) is strongly differentiable, has its range included in the domain of Hi_ l (uj) 
and converges strongly to one. Define B := Hi.l(uj) — A(a)Hx t L(u>)A(—a), or in 
other form: 

B = 2ta(p - cua) • Vp(- - x ) + a(Ap)(- - x ) - a 2 |Vp(- - x )| 2 . (4.72) 

A well known result says (see [S]) that || B[H\_l{u) + l] -1 / 2 || < const. 

By direct computation, i?(r) = B A{o)Wl{t, tu)A(— a); a rough estimate 
gives ||i?(r)|| < const(L)/y / r and even if the constant behaves badly with L, the 
norm is integrable in zero with respect to r. 

In conclusion: 

f P 

W L {(3,u)-S{f3) = - dTW L {l3-T,uj)R( T ) 
Jo 

P 

dTW L ((3-T,uj)BS(T), (4.73) 
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where the integral converges in norm. But uniformly in L and xo there exists a 
numerical constant such that: 

const 

\\W L {(3-T,u)B\\<a^L. (4.74) 
VP-r 

From l|4.74|l . 14.73JI and point i of Proposition El the needed estimate follows, 
ii. Using point i, the estimate l|4.3[) and the identity (a small enough): 

A(a)K ~ zW L {fi,u)]- l A{-a) = [f - zA(a)W L (/3, ^(-a)]- 1 = 

j>0 

■ {[A(a)W L ([3,u>)A(-a) - W L {J3,u)][t ~ zW L (f3, cu)}' 1 } 3 , 
the result follows. 

Corollary 1. The operator A(a)g UJ (£,, z; (3, t)A(— a) belongs to £?(L 2 ,L°°) if a 
is small enough, and uniformly in z, xn, u and L one has: 

i. \\A(a)g u (Z,z',P,T)A(-a)\\ < const(/3, t); 

ii. J A dye 2Q ly- x °l|r w (xo,y)| 2 < const(/3, r); 
hi. e Q l x ^l|T w (x,y)| < const(/3,T). 

Proof, i. It is an immediate consequence of Propositions 1 1 1 1 and 1121 

ii. Let <t> = A(a)g UJ (^, z\ [3,T)A(—a)ip, where 4> is bounded and continuous. 
From i it follows: 

|#xo)| = | / dy e Qp(y - Xo) 7Ux ,yMy)| < const(/3,r)||^|| L2 , (4.75) 

J A 

and the result follows from the representation theorem of linear and continuous 
functionals on L 2 . 

hi. Rewrite the identity <7^(£, z; j3, r) = <? w (£, z; (3, t/2)Wl(t/2, uu) in terms of 
integral kernels, use (|4.68|) . I|2.3|l . ii, the triangle and Schwartz inequalities, and 
the proof is completed. 

Let 5u> > be such that co = ojq + Soj 6 Q. Define the bounded operator 
g u (£, z; (3, t) given by the following integral kernel: 

T w (x,x') :=e^ (x < x,) T W0 (x,x'). (4.76) 

Equations (|4.3|l and (|4.58(1 imply that if 5uj is sufficiently small then there exists 
a numerical constant such that uniformly in L: 

supsup||[£- zS^/3)]- 1 !! < const (4.77) 

and: 

sup sup -zS* w {13)]- 1 - zWl^^)]-^ < const 5w. (4.78) 
We state now an important property of g u (£,, z; (3, t): 
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Proposition 13. Under the above conditions, there exists a numerical constant 
such that if Slu is small enough, then uniformly in £, z and L, the following 
B(L 2 ) estimate takes place: 

- zS u {P)]- x zS u {t) - /?,r)|| < const 8u. (4.79) 

Proof. The integral kernel of the operator [£ — zS w (0)]g w (^, z; 0, r) is given by: 

^(x.x')-* / dy 5 w (x,y;/3)f w (y,x'). (4.80) 

J A 

Let us notice a crucial property of the magnetic phase: 

p(x, y) + p(y, x') = p(x, x') + fl(x, y, x'), (4.81) 

where fl(x,y,x') =1/2B- [(y - x') A (x - y)]. Using and gSEI in 

we obtain: 

[£-zS„(/3)]s„(£,z;/3,t) = zS„(r) + J R, (4.82) 
where i? is an integral operator given by: 

_ z e ^(x,x') /" dy (e ^fl(x, y) x') _ 1} Ga)o i ( X) y ;/ J) Two ( y)X '). (4.83) 

Because |e J(5w fl ( x ^' x ) — 1| < Suj\x — y| |y — x'|, denoting with P the operator 
given by |x — y| |zG WOi i,(x, y; /3)| and with Q the operator corresponding to 

' y ~ *l[ r "°( y ' x ')|' }t follows that 1 1^1 1 ^ Ml^ll IIQH- Usin g E3>> CorollaryU] 
iii., H4.34fl and (|4.77|) . the proof is completed. 

Employing flOSJ) . (|4~79"1) . (|4~75|) and (|4~77|) in the next equality: 

g u (t,z; 0, 0/2) -g u (£,z; 0,0/2) = 
= [(£ - zW^M)- 1 - (£- 25 u (/3))- 1 ]^ L (/3/2,o;) + 
+ (£ - z^OS))- 1 *^^/^) - 5 w (/3/2)] + 

+ (£ - z5 w (^))- 1 ^( / 9/2) - 2; ft /3/2), (4.84) 

(|4.9(l is straightforward. 

Let us end this subsection by proving that the operator 

A iO (Z,z;0)=g u (t,z;!3,0/2)S u (J3/2) 

has the same trace as gu (£,z; 0, 0). Indeed, because A u fulfills the conditions 
of Proposition 1^1 and noticing that y(x, x') = — </j(x',x), one can write: 

trA,(£,z;/3) = f dxdx'T W0 (e,z;/3,/3/2;x,x')G W0 . L (x',x;/3/2) 

= / dxLT Uo (Z,z;0,0;xL,x)=tr 9uio (Z,z;0,0). (4.85) 
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4-.2. The proof of II and III. The analyticity of r oo {!3, z, u>o) in T> follows from 
the bound (see l|2.16[l ) |ff CT (C)| < const (a, K) \C\ where K is some compact in 
C\[l,oo). 

In what follows, we will prove that if z G T>q := {\z\ < 1}, then: 



lim r L (f3,z,w ) = roo(P,z,Uo). 

L — >-oo 



(4.86) 



Because \z\ < 1, the grand canonical pressure will be (see 1)2. H)^ : 

(4.87) 



Under the same conditions, the magnetization reads as 

~dW L 



r L <fi,z, Uo ) = — e J2-0 E s* 

n— 1 



du> 



(n/3,u ). 



(4.88) 



An important quantity is the integral kernel of the semigroup defined on the 
whole space: 



G Wi oo(x, x ; (3) — 



(2tt/3) 3 / 2 sinh(w/3/2) 



tanh (w/3/2) 



(e 3 A(x-x'))^ + (e 3 -(x-x')r 



(4.89) 
(4.90) 



Denote with: 



g(j3,w) = G aj , 00 (x,x;/3) 



)p/2 



(2tt/3) 3 / 2 sinh (w/3/2) " 



Then it is easy to see that if \z\ < 1: 



d In \ \^ z n g(n/3,uJo) r fa s e ^ z n dg 
Pco(P,z,uj q ) = — a . ^oo(P,z,u; ) = -— 2^ ~^( n ^' w o)- 



/3c ^ n duj 



One of the results in [M-M-P 1] can be adapted to our problem and gives: 



lim — tr 



dWr 



duo 



(4.91) 



If we prove the existence of a positive function / with at most polynomial growth 
such that: 

'dW L ' 



1 



duo 



(n(3,u; ) < f(nf3), 



(4.92) 



then 1)4.86)1 would be true. 

The next corollary is a direct consequence of Proposition [S]ii. 
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Corollary 2. Under the conditions of Proposition^ one has: 



lim -!- / dx [G Wo+ fc iL (x,x;/3) - G W0! l(x,x;/3)] 
5u\0 ouj J A 



(4.93) 



= tr 



dW L 



08, wo) = 2 / dT / rfxrfy 

f •/ A 2 



G W0 ,l(x, y;/3 - r)a(y - x) • (-iV y - w a(y))G Wo , L (y, x; r). 



Use g3BJ, and in and (jO^I follows. The proof of 

(|4.86|) is now completed. 

Remarks. 1. What can we say about the same problem for Fermi particles (say 
electrons, where Hi tL {u) should be replaced with the Pauli operator)? Knowing 
that in this case inf a(Hi t00 (uj)) = (i.e. is independent of w), the grand canoni- 
cal result (Lemma^) can be easily restated in terms of Fermi statistics: the only 
thing that changes is the domain on which the limit takes place i.e. C\ (-co, — 1], 
As for Theorem its proof was based on the fact that there exists a compact 
K C C \ [e' 3 "/ 2 , oo) which contains the circle centered in the origin with radius 
equal to Xl(/3, p,u>), for all L > Lq and for all strictly positive j3, to and p. For 
Fermi particles, if one fixes p and u> but makes (3 very large (lowers the temper- 
ature), then xl(P,p,oj) would be a very large positive quantity, therefore the 
above circle could intersect the negative cut. Of course, if the gas is diluted (say 
p and w fixed and (3 small), then it could happen that Xoo((3, p, w) < 1 which 
means that the circle never intersects the cut when L > Lq, therefore a similar 
proof can be provided. Our conclusion is that the extension of Theorem to a 
Fermi gas at low temperature is not trivial and remains an interesting problem. 

2. What about the higher derivatives with respect to u> (the susceptibility 
for example) at loq ^ 0? This also remains an open problem, even for the grand 
canonical ensemble. Nevertheless, we think that our approach (the modified per- 
turbation theory for Gibbs semigroups) could provide an answer to it. 
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